Exotic paired states with anisotropic spin-dependent Fermi surfaces 
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We propose a model for realizing exotic paired states in cold atomic Fermi gases. By using a spin dependent 
optical lattice it is possible to engineer spatially anisotropic Fermi surfaces for each hyperfine species, that are 
rotated 90 degrees with respect to one another. We consider a balanced population of the fermions with an 
attractive interaction. We explore the BCS mean field phase diagram as a function of the anisotropy, density, 
and interaction strength, and find the existence of an unusual paired superfluid state with coexisting pockets of 
momentum space with gapless unpaired carriers. This state is a relative of the Sarma or breached pair states 
in polarized mixtures, but in our case the Fermi gas is unpolarized. We also propose the possible existence of 
an exotic paired "Cooper-pair Bose-Metal" (CPBM) phase, which has a gap for single fermion excitations but 
gapless and uncondensed "Cooper pair" excitations residing on a "Bose-surface" in momentum space. 

PACS numbers: 74.20.-z, 74.25.Dw, 03.75.Lm 



Fermionic superfluidity is, with its many manifestations 
ranging from solid state physics to ultracold atom gases and 
neutron stars, one of the most fascinating phenomena. In re- 
cent years we have experienced a renewed interest in some of 
the basic questions behind this problem, motivated by techni- 
cal advances in the control and manipulation of cold atomic 
gases. [li] In particular, a great deal of attention is being de- 
voted to the quest for unconventional paired states of matter. 

Conventional superconductivity (BCS theory iQl) consists 
of pairing between fermions with opposite spin and equal but 
opposite momentum. In the presence of (partial) spin polar- 
ization more unusual superfluid states are possible. In the 
Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) statelH 0] pair- 
ing occurs across the mismatched Fermi surfaces so that the 
Cooper pairs have a finite center-of-mass momentum. |@] The 
excess unpaired fermions occupy a Fermi sea, so that there 
are gapless fermionic excitations coexisting with superfluid- 
ity. The FFLO state has been observed in solids only re- 
cently jslj. An interesting alternative to the FFLO state in par- 
tially polarized systems is astate with zero momentum pair- 
ing as suggested by Sarma With discon- 
nected regions of momentum space with pairing, separated 
by an unpaired polarized sea of fermions, the Sarma state is 
sometimes referred to as a "breached pair" phase ifioi [Till ). 
Recently a "pair density wave" or striped superconducting 
state with a spatially modulated order parameter similar to 
the FFLO state but without time-reversal symmetry breaking 
has been proposed to account for experimental observations 
in La2_^Ba:rCu04. 

These unconventional paired states might be accessible in 
cold atom experiments, with "polarized" Fermi gases in two 
hyperfine states having different populations. lH H 01 Ad- 
ditionally, one can explore the effects of differing masses in 
Fermi mixtures with two atomic species. lll6ll A mass "im- 
balance" between the two hyperfine states of a single atomic 
species can also be achieved in an optical lattice by tuning 
a spin-dependent hopping with light. Indeed, spin-dependent 



optical lattices have already been demonstrated lllTll . A proto- 
col to realize a Hubbard model with spin-dependent hopping 
has been proposed by Liu et al. in Ref. llSll . by tuning the lat- 
tice lasers between fine structure levels of ''"K atoms, and us- 
ing Feschbach resonances to manipulate the interactions. By 
varying the polarization and mass ratios it might be possible 
to realize unconventional paired states such as the FFLO or 
Sarma state. 

In this paper we suggest a different means to access un- 
conventional paired states in ultracold atomic systems. Our 
idea is to create mismatched Fermi surfaces in an unpolarized 
mixture. Consider an experiment with two hyperfine states 
(that we label t and [ hereafter) with equal population that are 
moving, for simplicity, on a two-dimensional square lattice. 
Imagine tuning the hopping so that one spin state hops pref- 
erentially along the .T-axis and the other preferentially along 
the j/-axis. We are interested in the situation where the Fermi 
surfaces of the two spin states are rotated by 90 degrees with 
respect to one another, and focus on a near neighbor hopping 
Hamiltonian with single particle dispersions, 

e^{k^,ky) = -2iaC0s(fc^) - 2ti,cos{ky) - /i, (1) 
ei{k,j:,ky) = -24 cos (fc^) - 2ta cos (ky) - /i, (2) 

with chemical potential /i. The ratio tt/ta < 1 determines the 
eccentricity of the two Fermi surfaces (see Fig[T]i. A similar 
Fermi surface geometry has been proposed for the ruthenate 
materials, | S^l arisin g in a multi-band model which pos- 
sesses spin nematic order 112 1 112211 . The nematic behavior orig- 
inates from hybridization of the bands and is a spontaneous 
symmetry breaking, while our Hamiltonian explicitly breaks 
time-reversal and rotational symmetry. 

As we discuss below, in the presence of an attractive inter- 
action this peculiar arrangement of anisotropic Fermi surfaces 
predicts a new class of gapless superfluid (G-SF) states within 
a BCS mean field theory, and might also give rise to an even 
more exotic paired but non-superfluid Bose-Metal phase. 

To be concrete we assume the particles interact through a 
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FIG. 1: (color online) Schematic Fermi surfaces for the disper- 
sions ecr{kx,ky) described by Eq.l|2j: (a) closed Fermi surfaces 
for small density; (b) open Fermi surfaces for higher density/large 
anisotropy. Arrows label two hyperfine states of the same atomic 
species. Dashed lines delimit the pockets of the momentum distribu- 
tion in the gapless G-SF solution. 



short-range s-wave potential, that we represent using the at- 
tractive Hubbard model: 



H 



k.a 



Ck., 



(3) 



where cj. ^ creates a fermion with spin cr = f, j at momentum 

k, Tiiu = cl^Cia is the local on-site density, and U is the inter- 
action strength that we will take to be negative (attractive). 

By simple inspection of the spin-dependent Fermi surfaces 
in FiglT] we can anticipate several different pairing possibil- 
ities within a BCS treatment. For a very strong attractive in- 
teraction, \U\ >> ta,ti„ a state with zero momentum pairing 
and a fully gapped Fermi surface is expected. For smaller U a 
zero-momentum paired state with gapless single fermion ex- 
citations analogous to the Sarma state in polarized mixtures 
could occur. Alternatively, pairing could occur across the two 
mismatched Fermi surfaces with a finite momentum conden- 
sate, an unpolarized analog of the FFLO state. 

We now implement a self-consistent BCS approximation 
to obtain the mean-field phase diagram and explore the prop- 
erties of the resulting paired superfluid states. A discussion 
of more exotic non-superfluid paired phases which are inac- 
cessible within mean field theory will be presented below. 
We consider an on-site pairing of the form -~U/V{cjiCj^) = 
A exp (iQj) with gap parameter A. Here V — is the num- 
ber of sites in the system and j denotes the center of mass po- 
sition of a pair in real space. The resulting quadratic Hamilto- 
nian can be readily diagonalized, yielding quasi-particle exci- 
tations with dispersion: 



(4) 



where we have defined, 

e±(fc) = h(-fc+|)±q(fc + |)]/2. (5) 

Consider first the case of zero pairing momentum Q = 0. 
Besides the normal phase with A = 0, we expect a fully 



gapped BCS-like solution when E±{k) > for all momenta 
in the Brillouin zone (BCS state). This occurs when the at- 
tractive interaction is sufficiently strong. For smaller A a gap- 
less superfluid phase (G-SF) is also possible. [l^ Tsl Gapless 
fermion excitations occur at momenta that satisfy E^ik) ~ 
or E-{k) = 0, or equivalently. 



eT(fc)ex(fc) = -A2 



(6) 



The resulting closed curves in momentum space are sketched 
schematically as dashed lines in Fig. 1 . In the case of closed 
Fermi surfaces as shown in Fig. 1(a) there are four small pock- 
ets near the Fermi points (ifcj^^, 0) and (0, ±fc]^^). The so- 
lution with open Fermi surfaces (see Fig. 1(b)) has two large 
pockets that wrap around the Brillouin zone. 

It is instructive to compute the momentum distribution 

Ck^a)- There 



function for the two spin species, na{k) 
are three distinct regions in momentum space determined by 
the conditions: (i) E^{k) > and E+{k) < which occurs 
inside the f pockets; (ii) E+{k) > and E^{k) < which 
occurs inside the j pockets; or (iii) both E± (fc) > which oc- 
curs over the rest of the Brillouin zone. For momentum where 
(i) is satisfied we obtain n|(fc) = = 0, whereas in 

region (ii) one has n^{k) = 0,nj^(fc) = 1. Thus, inside these 
pockets the fermions are completely unpaired, and there is no 
mixing between the two spin species. Finally, in region (iii) 
the fermions are completely paired and one obtains; 




6+(fc) 



Ve+(fc)2 + A2 



(7) 



The equivalence of the two momentum distribution functions 
and the restoration of the four fold symmetry of the square 
lattice are presumably artifacts of BCS mean-field theory, and 
would be modified when fluctuations are taken into account. 
Even within BCS mean-field theory, the spectral functions of 
the two fermion species are different from one another 

While possessing regions in momentum space with com- 
plete pairing and other regions with no pairing whatsoever is 
reminiscent of the Sarma phase, this novel unpolarized gap- 
less superfluid (G-SF) phase has unpaired fermion pockets of 
both spin species. 

We now analyze the (mean-field) stability of the normal, 
BCS and gapless superfluid (G-SF) phases, by solving the 
BCS saddle point (gap) equations numerically on the square 
lattice. We work in the canonical ensemble, with fixed density, 
which is the relevant situation for cold fermionic gases, also 
allowing for finite momentum pairing with momenta Q ei- 
ther along the x, y-axes or along the diagonals in the Brillouin 
zone. The Helmholtz free energy in the canonical ensemble 
is computed and minimized in systems sizes up to 300 x 300 
sites following the prescription detailed in Ref.|23;]. 

For small values of the attraction \ U\/twe find that the sys- 
tem is in the normal state, independently of the total fermion 
density, n, or the value of the anisotropy tb/ta, while for 
large values of the interaction strength, the system is in a fully 
gapped superfluid state with Q = (BCS). The interesting 




FIG. 2: Mean-field phase diagram of the attractive Hubbard model as 
a function of the anisotropy tt jta and fermion density n, for values 
of the interaction (a) U — — 3.5t and (b) U = —4t. We find regions 
occupied by a fully paired solution with Q — (BCS), normal (N), 
gapless superfluid (G-SF), and phase separation (PS). 



regime is for values of 2t < \U\ < At, where we find a com- 
petition of phases. In Fig|2]we show the phase diagram for 
two different values of \U\/t, as a function of the fermion 
density and the anisotropy ti,/ta- At \U\/t = 3.5, the normal 
phase occupies a region at high densities and high anisotropy, 
while the superconducting phase (BCS) dominates at inter- 
mediate and small anisotropy. Separating these two phases 
is a window where the G-SF state is realized. There is also 
a tiny sliver where phase separation takes place, that coin- 
cides with the onset of the instability toward the G-SF state. 
We determined the boundaries of this region by identifying 
densities with negative compressibility dn/d^i, and using the 
Maxwell construction. Contrary to the Sarma phase in po- 
larized mixtures, the G-SF state is stable in a wide region of 
parameters . 1 1 2[ 1 1 311 At larger values of the interaction the nor- 
mal phase becomes energetically unfavorable, leaving a re- 
gion of the G-SF phase that shrinks with increasing interac- 
tion. In all cases the G-SF phase only occurred in the with 
open Fermi surfaces/large pockets (see Fig. 1(b)). We also 
looked for energy minima at finite momentum, but we have 
not found solutions for the system sizes and in the parameter 
space that we studied. 

We now turn to the possibility of paired states which are 
not superfluids. This is motivated by recent work|j3 which 
introduced and studied a hard core boson model with ring 
exchange, and made strong arguments for the existence of a 
novel Bose-Metal phase, a D-wave Rose Liquid. Here we 
will refer to this phase as a D-wave Bose-Metal (DBM). The 
proposed Bose-Metal phase is not a superfluid, having zero 
condensate and superfluid densities, but possesses gapless ex- 
citations which reside on surfaces in momentum space. In 
contrast to a conventional Fermi liquid, in the DBM phase the 
gapless excitations cannot be described in terms of weakly in- 
teracting quasiparticles. The low energy excitations are intrin- 
sically strongly interacting fluctuations, and a single particle 
description is incorrect. In the DBM phase pairs of bosons 
possess D-wave correlations, but the phase is not a paired- 
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FIG. 3: Processes that would give origin to an effective ring- 
exchange of pairs, fermions of different hyperfine species (with op- 
posite spin) prefer to move along perpendicular directions (left). 



boson condensate. 

A connection between our model of spin-fuU fermions with 
spin-dependent spatially anisotropic Fermi surfaces and the 
Bose-Metal phase explored in Ref . t24i1 can be made as fol- 
lows. First we define an on-site "Cooper pair" operator. 



h] = 



This operator is a "hard-core" boson, commut- 
ing at different sites, = for i ^ j, yet satisfying 

blbl = 0. Now consider the regime with \ U\ >> ta,ti, where 
all of the fermions are tightly bound into on-site Cooper pairs. 
The goal is to derive an effective Hamiltonian for these hard- 
core bosons, by considering a perturbation expansion in pow- 
ers of ta/|J7|, it/ 1?7|. At second order one obtains a boson 
hopping term Hj with strength J tatb/\U\, 



H.I 



ihj) 



H.c. 



(8) 



The hoping J is proportional to a product of ta and tf, since 
it is necessary to hop both the t-spin and J,-spin fermion in 
order to move the hard-core boson. Notice that in the extreme 
anisotropic limit with ti, ~ and 7^ it is actually impos- 
sible for a single boson to move, J = 0. 

But consider the boson ring-exchange term as depicted in 
Figl3] which involves two bosons on opposite corners of an 
elementary square plaquette rotating by ±90 degrees. 



i?ring = K^ing ^ ^1 62^3^4 + H.C, 



(9) 



plaqucttcs 



with i = 1,2,3,4 labeling sites taken clockwise around 
the square plaquette. The ring-exchange coupling strength 
-K^ring > can be computed at fourth order in ta.b/\U\. In 
the highly anisotropic limit with >> it, the process de- 
picted in Figl3]dominates, giving K^i^g ^ t'^/\U\'^. Remark- 
ably, while the single boson hopping term vanishes in the 
extreme anisotropic limit, tb/ta — » 0, the ring term which 
hops pairs of bosons is non-zero. Thus, with increasing 
anisotropy, ta/tb, the ratio /^ring/ J ~ {ta/tb)^ increases, and 
the ring term becomes increasingly important in the Hamilto- 
nian, Hb = Hj + i/ring- 

Consider the phase diagram for the boson Hamiltonian Hb, 
when the bosons are at some generic density, Ub = {b\bi). For 
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^ring « J one expects a Bose-condensed superfluid phase. 



In Ref. II24I1 variational wavefunction studies found that the D- 
wave Bose-Metal (DBM) was stabilized for A',ing/ J > 2 and 
for boson densities larger than ji^ > 0.4 (or, under particle- 
hole symmetry, smaller than 0.6). At other densities in this 
large i^ring regime, phase separation into a (probable CDW) 
phase at half-filling and a zero density state was found. 

The boson superfluid state corresponds to a fully gapped 
(conventional) BCS state of the fermions. However, the DBM 
phase of the bosons would correspond to a truly exotic non- 
superfluid paired state of fermions, with a full gap for the sin- 
gle fermion excitations, but with uncondensed "Cooper pairs". 
The Cooper-pairs would form a collective state with gapless 
excitations along surfaces in momentum space - a "Cooper- 
pair Bose-Metal" (CPBM)! 

Since the fermion density, n = n| + = 2nb, this would 
suggest that the CPBM phase might be present in the phase 
diagram in the regime of intermediate U /t, strong anisotropy 
with tb/ta small and for densities near half-filling, n > 0.8. 
As one can see from Fig. 3(b), this is precisely the regime 
where the BCS mean-field treatment finds the stable gapless 
superfluid (G-SF) state. The G-SF and CPBM phases are dra- 
matically different. The G-SF state is a superfluid with gapless 
fermion excitations, while in the CPBM phase the fermions 
are fully gapped but the state is not a superfluid. Which of 
these two states, if either or perhaps both, is present in the 
true phase diagram of the original attractive Hubbard model 
is unknown, and certainly worthy of future investigation. Re- 
cent studies of the boson ring model on the two-leg ladder 
have found compelling evidence for a quasi- ID Bose-Metal, 
a ladder descendant of the 2D DBM state. Q 

To summarize, we have proposed a model that allows for re- 
alizing exotic paired phases of unpolarized fermion mixtures 
using spin-dependent optical lattices. The main ingredient is 
the existence of mismatched Fermi surfaces between both hy- 
perfine states. We explored the BCS mean-field phase dia- 
gram for two anisotropic Fermi surfaces rotated by 90 degrees 
respect to one another, checking for the relative stability of 
the BCS, FFLO, gapless superfluid and normal metal phases. 
Other mean-field states such as a nematic superfluid l'26'l, or 



multi-modal FFLO states. 11271 12811 could be considered in fu- 
ture work. We also argued that anisotropic Fermi surfaces plus 
attractive interactions leads to an effective model of Cooper- 
pairs with a ring-exchange term, that may allow to realize 
a paired but non-superfluid Bose-metal phase. Whether this 
Cooper-pair Bose-Metal phase appears for the attractive Hub- 
bard model will be explored in future work. The Fermi surface 
geometries presented here are just one example that demon- 
strate the rich physics that emerges once one allows for spin 
dependent Fermi surfaces. Other possibilities can also be ex- 
plored, such as a combination of circular and elliptical Fermi 
surfaces, with the only requirement that both Fermi surfaces 
enclose the same area. 
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